Damped Bloch Oscillations of Bose-Einstein 
Condensates in Disordered Potential Gradients 



S. Drenkelforth\ G. Kleine Biiningi, J. Will^, T. Schulte^ 

N. Murray^ W. Ertmer^ L. Santos^, and J.J. Arlt^ 

^ Institut fiir Quantcnoptik. Leibniz Univcrsitat Hannover, Welfengartcn 1, 
D-30167 Hannover. Germany 

^ Institut fiir Theoretische Physik, Leibniz Universitat Hannover, Appclstrafie 2, 
D-30167 Hannover, Germany 

E-mail: drenkelf orthSiqo . uni-hannover . de 

Abstract. We investigate both experimentally and theoretically disorder induced 
damping of Bloch oscillations of Bose-Einstein condensates in optical lattices. The 
spatially inhomogeneous force responsible for the damping is realised by a combination 
of a disordered optical and a magnetic gradient potential. We show that the 
inhomogencity of this force results in a broadening of the quasimomentum spectrum, 
which in turn causes damping of the centre-of-mass oscillation. We quantitatively 
compare the obtained damping rates to the simulations using the Gross-Pitaevskii 
equation. Our results are relevant for high precision experiments on very small forces, 
which require the observation of a large number of oscillation cycles. 
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1. Introduction 

The ability to realise ultracold quantum gases in periodic and disordered potentials has 
enabled detailed studies of fascinating effects originating in solid state physics. Ongoing 
investigations of single particle phenomena such as Anderson localisation as well as many 
particle effects like the Bose-Glass phase and the Mott insulator [H Ej [3], HI [5l [6] show 
the variety of possibilities these ensembles offer. Especially the non-intuitive dynamics 
of quantum gases in periodic potentials is of interest for theoretical and experimental 
investigations, since quantum gases have enabled the first direct observation of Bloch 
oscillations [8] in tilted periodic potentials P, [ini [H] . In these systems the periodicity 
leads to an oscillatory motion instead of a linear acceleration of the particles subjected 
to an external force. 

In solid state systems scattering at imperfections of the crystal structure leads to 
damping of Bloch oscillations on timescales much shorter than the oscillation period 
itself. Therefore Bloch oscillations of electrons are only observable in semiconductor 
super lattices [12], where the large spatial period leads to high oscillation frequencies, 
which are faster than the damping. Optical lattices on the other hand constitute perfect 
optical crystals and allow for the observation of long lived Bloch oscillations [131 El CSl 
[T6] . The experimental control of lattice parameters such as lattice depth and spacing, 
the possibility to detect the atomic cloud with absorption imaging and the very small 
momentum spread of Bose-Einstein condensates (BEG) have enabled detailed studies 
of this quantum effect. 

A comparison of these systems gives rise to the question how the controlled addition 
of disorder to an optical lattice will affect the dynamics of particles in such a periodic 
potential. Disorder can be realised with additional optical potentials [HI [181 [IS El EDI 
EH ESI El [7], impurity atoms [231 El] or the roughness of the trapping potential close to 
the surface of atom chips [25]. The simultaneous application of a homogeneous force and 
the optical disorder potential constitutes a spatially inhomogeneous acceleration. This 
inhomogeneity can have important consequences for the application of Bloch oscillations 
as a sensitive tool for high precision measurements of small forces [151 EHl EH El [13] . 
since it leads to a dephasing of the quasimomentum and thus to a damping of the centre- 
of-mass oscillation [28]. Therefore a detailed quantitative understanding of the effect of 
the disorder and the underlying mechanism is indispensable for future applications. 

We investigate the effect of a small disordered potential on Bloch oscillations of 
Bose-Einstein condensates in a ID optical lattice potential. Our results indeed show 
that the inhomogeneity leads to significant damping of the centre-of-mass motion. The 
damping rate increases with disorder depth and is quantitatively compared to numerical 
simulations using the Gross-Pitaevskii equation (GPE). Furthermore we show that the 
disorder induced broadening of the quasimomentum distribution, which is the underlying 
mechanism for the damping of Bloch oscillations, reduces the fraction of atoms in the 
BEG. 
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2. Bloch oscillations in periodic potentials 

The acceleration of particles in periodic potentials leads to an oscillatory motion instead 
of a linear increase in velocity. Since this is a pure single particle effect, it is possible 
to describe the underlying physics in a ID model, while a quantitative analysis of the 
effects of disorder and interactions require a 3D description, as reported in a previous 
work [28] . We briefly discuss the main features necessary for the understanding of Bloch 
oscillations while a comprehensive review can be found in [29] . 

The periodicity of the potential implies that the eigenf unctions obey the same 
translational symmetry as the potential 

<P{z + d)=e'^U{z), (1) 

where the phase difference q from site-to-site is called quasimomentum and d is the 
lattice constant. Since the eigenf unctions are periodic, it is possible to restrict the 
description of the dynamics to the first Brillouin zone According to the 

acceleration theorem q{t) = g(0) + ^ t an additional potential gradient will cause 
the quasimomentum to evolve linearly in time, because the energy offset from site- 
to-site leads to a linear increase of the phase difference from site-to-site during the time 
evolution. In combination with the periodicity of the band structure this causes an 
oscillatory motion, since the group velocity is proportional to the derivative of the band 
structure. The oscillation period T = 2'jTh/F d depends on the applied force F and the 
lattice constant. The resulting amplitude zbo = ^/2|-F| is given by the width of the 
first band A and the force. Since this amplitude is only a few micrometres for typical 
experimental parameters, the Bloch oscillations are conveniently analysed in momentum 
space via time-of-fiight (TOF) absorption imaging. 

While the oscillation can be described in a simple ID model, the complex dynamics 
in a disordered potential gradient and the role of interactions have to be analysed with 
a full 3D model. All numerical simulations presented in this work are obtained using 
the 3D Gross-Pitaevskii equation 



*2 

-— + Vl{z) + VMF{r) + Vgradiz) + 9 |^(r)|' 



*(r) = /i^(r). (2) 



The cylindrically symmetric magnetic trapping potential is VMpiT^) = hmUpp'^ 



with trapping frequencies Up and oo^- The optical potential is Vl{z) 



1 2 2 

s Er cos^ (kz), where s is the lattice depth in units of the recoil energy 

and /c = 27r/A is the wave vector of the optical lattice. The additional homogeneous 

gradient potential is given by Vgrad{z) = F z. 

Due to interactions between the atoms, Bloch oscillations of Bose-Einstein 
condensates suffer, contrary to the case of fermions [13], from dynamical instabilities [311 
[32] . In the outer half of the Brillouin zone the nonlinear coupling leads to an exponential 
growth of small perturbations. Independent of inhomogeneities in the potential they 
are responsible for a damping of the Bloch oscillations. Therefore it is necessary to 
reduce the dynamical instability in order to investigate the effect of the disorder on the 
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Figure 1. Absorption images after a TOF of SOms, for varying Bloch oscillation times. 
The lattice depth was 2Er and the acceleration 2.4m/s^. 



Bloch oscillations. This can be accomplished with a combination of a high potential 
gradient [35] and a reduction of the nonlinear interactions. In recent work long lived 
Bloch oscillations of Bose-Einstein condensates were realised p!6l [Hj by decreasing 
the s-wave scattering length with a Feshbach resonance, essentially turning off the 
atomic interaction. In the experiments reported here the interactions were reduced 
by decreasing the density of the BEC, to enable the observations of disorder induced 
damping of Bloch oscillations. 
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Figure 2. Ccntrc-of-mass position of Bosc-Einstcin condensates after a TOF of 30ms, 
as a function of the Bloch oscihation time. The parameters were identical to Fig. [TJ 



3. Experimental realisation 

The experiments were performed with ^^Rb Bose-Einstein condensates in the F = 
2, mp = 2 state. A detailed description of our apparatus is given in [33]. Nearly- 
pure BEC of up to = 3 X 10^ Atoms are produced, however all experiments described 
here were carried out with A^ = 5 x 10^ atoms to reduce the interaction energy, and to 
ensure that no significant thermal background is present. After production of the BEC 
the magnetic offset field was adiabatically increased to 91G within 440ms, thus lowering 
the radial trapping frequency to u± = 271 ■ 29Hz. This reduces the interaction energy 
by a factor of 6.25, which in combination with low atom numbers sufficiently inhibits 
the dynamical instability, such that up to four undamped Bloch oscillation periods 
can be observed. The following procedure was used to observe the Bloch oscillations. 
After decreasing the radial trapping frequency the intensity of the optical lattice was 
adiabatically increased to its final value within 60ms. Subsequently the atoms were 
subjected to either a homogeneous potential gradient or the spatially inhomogeneous 
one for a variable time. Finally all potentials were turned off at the same time and the 
atomic cloud was detected after a time-of-fiight of 30ms by absorption imaging. 

The ID optical lattice is provided by a standing light field at a wavelength of 
825nm, which is superimposed on the axial direction of the magnetic trap with a waist 
of c^o = 140/im at the position of the atoms. The investigations were performed at a 
lattice depth of 2Er, since low lattice depths lead to a large width of the energy band 
and therefore to a high maximal group velocity. This results in an oscillation amplitude 
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Figure 3. Centrc-of-mass oscillation position of a BEC for a lattice depth oi 2 Er, an 
acceleration of 2.4ni/s^ and disorder depths of (black, solid line). 70 x 10~^ E^ 
(red, dotted line) and 130 x 10~^ Er (blue, dashed line), obtained from a numerical 
solution of the GPE. 

of the centre-of-mass motion of 80yU.m after a TOF of 30ms, which can easily be detected, 
while the Bloch oscillations are not affected by Landau-Zener tunnelling [M] . 

The homogeneous potential gradient is provided by magnetic coils in anti-Helmholtz 
configuration, which produce gradients of up to 3.7G/cm. This corresponds to an 
acceleration of 2.4m/s^. 

The inhomogeneity is realised by a disordered optical dipole potential, generated by 
imaging a randomly structured chrome substrate radially onto the BEC, as described 
in a previous publication [2]. The correlation length of the disorder is 8yum and its 
depth was varied between and 135 x 10"^ Er, where the depth is defined as twice the 
standard deviation analogue to [T8] . 

Figure [1] shows absorption images of Bose- Einstein condensates undergoing Bloch 
oscillations after 30ms TOF for a lattice depth of 2Er without disorder. One clearly 
recognises the oscillation as a motion of the central peak and the periodic appearance 
of a second peak replacing the main one during an oscillation cycle. Figure [2] shows the 
centre-of-mass oscillation as well as a theoretical prediction based on a band structure 
calculation without free parameters. The position of the Bose- Einstein condensate zbo 
was calculated by numerically summing the weighted axial positions zbo = J2z w^^- 

total 

The measured Bloch period was 4.5ms and the oscillation amplitude was 73/xm which 
is in good agreement with the theoretical calculations of 4.66ms and 80/im for an 
acceleration of 2.4m/s^ and a lattice depth of 2Er. 
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Figure 4. Ccntrc-of-mass position of Bose-Einstein condensates performing damped 
Blocli oscillations. The disorder depths were 35, 105 and 135 x 10"^-^^; from top 
down. A fit to the data (solid red line) and the Gaussian envelope (dashed black line) 
due to the damping are also shown. Note that the recorded oscillation time is reduced 
for increased disorder depth since the broadening of the quasimomentum spectrum 
causes a strong reduction of the contrast (see Fig. [8|). This significantly reduces the 
signal-to-noise ratio of the absorption images. 

The theory of disorder induced damping of Bloch oscillations was investigated in 
detail in a previous publication [28] and we only briefly review the important features 
for the interpretation of the experimental results. 

For the analysis of Bloch oscillations in an inhomogeneous potential gradient the 
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Figure 5. Damping coefficient 1/cr of the centre-of-mass oscillation in the disordered 
lattice potential. The damping coefficients were obtained by applying the same 
fit procedure to the experimental data of Fig. 2] and the results of the numerical 
simulations. The red dots represent the experimental data and the black dots the 
simulations. Similar to the experimental case, the number of Bloch oscillation periods 
used in the fits to the simulations was reduced for increasing disorder depth. The 
shaded area corresponds to atom numbers used in the simulations ranging from 3.5x10^ 
to 6.5 X lO'^. 



GPE (equation [2]) has to be modified by expanding the acceleration term to 

Vgrad = F Z + Vdis{z), (3) 

where Vdis{z) denotes the additional optical potential which constitutes the disorder. 

Figure [3] shows the centre-of-mass position obtained from numerical simulations of 
damped Bloch oscillations for a typical disordered potential used in the experiment. One 
clearly recognises that the damping of the oscillations strongly depends on the disorder 
depth. The damping can be understood qualitatively in terms of the evolution of the 
phase difference from site-to-site. For an undisturbed homogeneous potential gradient 
it develops in time according to 

6E 

A0(t) = — t, (4) 

where SE is the energy offset and A(f) the time dependent phase difference between 
neighbouring sites. 

The disorder gives rise to a spatially varying energy difference from site-to-site 
SE{z). Therefore the phase evolution and the quasimomentum vary across the lattice. 
This broadening of the quasimomentum spectrum causes the damping of the Bloch 
oscillations, since each q corresponds to one group velocity. 
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Figure 6. Grey scale plot of the evolution of the momentum distribution during a 
Bloch oscillation for a disorder depth of 105 x 10^^ Er- The lattice depth and the 
acceleration were identical to the case of Fig. [3l 
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Figure 7. Absorption images of Bose-Einstein condensates for an oscillation time of 
0.5 Bloch periods at an acceleration of 2.4m/s^ and a lattice depth of 2 Er after a 
TOF of 30ms. The left column shows the case without added disorder, whereas the 
right column corresponds to a disorder depth of 105 x 10""^ Er- The top row shows the 
absorption images, the middle row corresponds to the associated axial density profiles 
and the bottom row contains the momentum spectrum obtained from our simulation. 



Figure H] shows the experimental observation of disorder induced damping of Bloch 
oscillations for various disorder depths. Note that the lattice depth and the acceleration 
are identical to the undamped oscillation in Fig. O The graph clearly shows the distinct 
reduction of the oscillation amplitude for increased disorder. The solid lines are a fit 
to the data, with the damping coefficient and the periodicity as free parameters. We 
generate the fit function for the damped Bloch oscillation znBoif) by multiplying the 
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Figure 8. Fraction of atoms in the BEC after a time-of- flight of 30ms as a function 
of the Bloch oscillation time. Shown are four depths of the disorder potential of QEr 
(black circles), 35 x IQ-^Er (red squares), 70 x lO'^Er (blue dots), and 135 x 10"^^;^ 
(green triangles). The curves correspond to a fit with an exponential decay. Residual 
imaging effects (see figure [7]) result in an overestimate of the total number of atoms. 
Therefore Njrac is 0.5 even for pure BEC without any discernible thermal fraction. 

time evolution of the undamped oscillation ZBoi^)-, obtained from the band structure 
calculation, with a Gaussian envelope 

ZDBoit) = Ao e-'^l"^ ZBoit). (5) 

The shape of the envelope is a consequence of the broadening of the quasimomentum 
distribution in time. The width of this distribution determines the oscillation amplitude, 
which is reduced compared to the single particle picture, because separate parts of 
the ensemble simultaneously undergo different phases of the oscillation. Based on the 
assumption that the width of this distribution increases linearly in time, a Gaussian 
envelope of the damping amplitude is expected [SHI ED] • 

To quantitatively compare the experimental data to the numerical solutions of the 
GPE, we show the resulting damping coefficients as a function of the disorder depth 
in Fig. [5l The parameters of the simulations correspond to the experimental ones 
for a typical realisation of the disordered potential. Since the damping rate strongly 
depends on the depth of the disorder potential, the horizontal error bars represent 
an experimental uncertainty in this depth of 25%. This was estimated by evaluating 
the depth of the used disorder potential at different positions, while small deviations 
between the exact shape of the disorder potential in the experiment and in the numerical 
simulation were not accounted for. The shaded area corresponds to an uncertainty in 
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Figure 9. The decay times r of Fig. [8] are shown as a function of the disorder depth. 
Strong reduction is observed even for moderate disorder depths. 



the atom number of 30%, accounted for in the simulations. Within this uncertainty, we 
observe good agreement between experimental and theoretical values of the damping 
coefficients. 

To show, that the broadening of the quasimomentum spectrum is the underlying 
mechanism of the damping, we analyse the experimental data with a second approach. 
Since the dephasing of the quasimomenta increases the number of atoms with varying 
momenta, the number of atoms in the BEC peaks is reduced [35]. This picture is 
confirmed by Fig. [6l which shows numerical results for the time evolution of the 
momentum spectrum in the presence of a disordered potential gradient. The disorder 
induced dephasing leads to an increased blurring of the sharp momentum distribution, 
as the quasimomenta undergo a non-uniform evolution. 

This behaviour is consistent with the absorption images shown in Fig. [71 In the 
undisturbed case sharp BEC peaks are visible whereas the disordered case exhibits 
a clearly discernible background and a broadening of the peaks due to the dephased 
quasimomenta. 

To analyse this effect quantitatively, we estimate the number of atoms in the BEC 
iVBEC at different times of the Bloch oscillation by fitting Thomas-Fermi profiles to 
the characteristic peaks of the BEC. The fraction of atoms in the BEC is calculated 
by comparing the number of atoms in these peaks with the total atom number 

Nfrac = NBEc/Ntotal- 

Figure [HI shows this fraction as a function of the Bloch oscillation time for various 
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depths of the disorder potential. The hnes are fits to the data with an exponential decay. 
A clear reduction in the occupation of the BEC peaks during the evolution of the Bloch 
oscillation is visible. In all cases the dephasing of the quasimomenta precedes the onset 
of the damping shown in Fig. HI This is in agreement with a Gaussian shape of the 
damping envelope given in Eqn. [5]and confirms that the dephasing of the quasimomenta 
is the underlying mechanism for the damping of the centre-of-mass motion. Hence only 
a significant broadening of the quasimomentum leads to a relevant damping. 

Note that even at a disorder depth of OEr a reduction and therefore a dephasing 
is observed due to the interparticle interactions, while the centre-of-mass oscillation in 
Fig. mis still unaffected. 

The decay times r from Fig.[8]are shown in Fig.[9]as a function of the disorder depth. 
The decreasing fraction of atoms in the BEC peaks again confirms that stronger disorder 
leads to faster dephasing of the momentum distribution and that this broadening of the 
quasimomentum is the underlying mechanism for the damping of the Bloch oscillations. 

5. Conclusion 

We have presented the first experimental investigation on disorder induced damping 
of Bloch oscillations of Bose-Einstein condensates. The application of an additional 
disorder potential during the oscillation leads to a strong damping of the centre-of-mass 
motion and to a significant reduction of the fraction of atoms in the BEC. The observed 
damping rates are in good agreement with predictions based on numerical solutions of 
the full Gross-Pitaevskii equation and show that the underlying physical mechanism for 
the damping is the broadening of the quasimomentum spectrum due to the spatially 
varying phase evolution of the condensate. 

We show that even very small disorder results in fast dephasing of the 
quasimomentum and therefore damping of the Bloch oscillation. Since the disorder 
presented here is equivalent to a spatially inhomogeneous force, the results are of special 
interest for the application of Bloch oscillations for high precision spectroscopy of very 
small forces. To reach high precision in such experiments it is essential to follow a large 
number of Bloch oscillations. This number may be reduced if the observed force is 
spatially inhomogeneous on length scales comparable to the extend of the condensate. 
The good agreement between theory and experiment shows the applicability of our 
method to analyse the effects of spatially varying forces, and allows for estimates of the 
effect of small inhomogeneities for future experiments. 
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